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The average number Srf(t) of distinct sites visited up to time t by N noninteracting random 
walkers all starting from the same origin in a disordered fractal is considered. This quantity SN{t) 
is the result of a double average: an average over random walks on a given lattice followed by an 
average over different realizations of the lattice. We show for two-dimensional percolation clusters 
at criticality (and conjecture for other stochastic fractals) that the distribution of the survival 
probability over these realizations is very broad in Euclidean space but very narrow in the chemical or 
topological space. This allows us to adapt the formalism developed for Euclidean and deterministic 
fractal lattices to the chemical language, and an asymptotic series for Sjv(i) analogous to that 
found for the non-disordered media is proposed here. The main term is equal to the number of 
sites (volume) inside a "hypersphere" in the chemical space of radius L[ln(N) /c] 1 ^ where L is the 
root-mean-square chemical displacement of a single random walker, and v and c determine how 
fast 1 — Tt(l) (the probability that a given site at chemical distance I from the origin is visited 
by a single random walker by time t) decays for large values of l/L: 1 — Tt(() ~ exp[—c(£/L) v ]. 
The parameters appearing in the first two asymptotic terms of S'jv(t) are estimated by numerical 
simulation for the two-dimensional percolation cluster at criticality. The corresponding theoretical 
predictions are compared with simulation data, and the agreement is found to be very good. 

PACS numbers: 05.40.Fb,05.60.Cd,66.30.Dn 



I. INTRODUCTION 

In the preceding paper |IJ we studied the problem of evaluating the average number S^it) of distinct sites (or 
territory covered) by N independent random walkers that, all starting from the same site, diffuse on a deterministic 
fractal during time t. An answer to this problem was given in Ref. M in terms of an asymptotic series completely 
analogous to that given in Refs. |^|| for Euclidean substrates. In this present paper we turn our attention to the 
case of stochastic (or disordered) fractal media and again propose the same kind of asymptotic solution, although, in 
contrast with all the previous cases, we now have to translate the procedure to the topological or chemical language. 

Stochastic media are not constructed by the iteration of an invariable rule, such as that used in deterministic fractals, 
but are rather the result of a random process. Consequently, the concept of generator and rigorous self-similarity 
is absent and their fractal nature is recognized by the scaling of statistical quantities. Many natural objects share 
this statistical-fractal structure so that stochastic models seem to be more suitable to represent diffusion in real 
media. Of particular interest is the percolation model that has been used to characterize many disordered systems 
. This model is constructed by filling a regular lattice with "occupied" sites with a certain probability. Nearest 
neighbor occupied sites are supposed to be connected and form a series of clusters. At a certain critical concentration 
p c an infinite cluster appears, which is called the incipient percolation aggregate or percolation cluster at criticality. 

The quantity we are interested in, 5jv(t), is, for disordered media, the result of a double average: an average over 
the walks that the N random walkers can perform over a given lattice, followed by an average over many (ideally, 
all) realizations of the random lattice. This situation leads to certain subtleties, absent in deterministic fractals, that 
require special treatment. In particular, 5jv(i) can be expressed by 

M*) = (E {i-[r t (r)f }) , (1) 

where the sum is over all the sites of each fractal lattice, T t (r) is the survival probability, i.e., the probability that site 

r has not been visited by time t by a single random walker starting from the origin, Y] (l — [r t (r)] Ar | represents the 

mean territory explored by the N random walkers on a given lattice (the first average) ||, and ([• • •]) indicates that 
the average (the second average) of [■ • •] has to be performed over all possible stochastic lattices compatible with the 
random generation rules. Equation (|l|) can be rewritten as 
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oo ln(m) \ 

^(t) = E(E {l-lUrm,i)f}) , (2) 

m=0 \ i=l / 

where r m j stands for the i-th site out of n(m) that are separated from the origin by a Euclidean distance between 
mAr = r m and (to + l)Ar with Ar small (say, of the order of the lattice spacing). If r t (r mj i) is almost independent 
of i and the lattice realization, i.e., if the fluctuations in the probability density r t (r m i ) follow a narrow distribution, 
then one could approximate r t (r m i ) ~ (T t (r m i )} = T t (r m ), and therefore estimate S]y(t) by 

oo 

S N (t)=Y,{l-Ft(r m )} N }(n(m)) , (3) 

rn— 

where (n(m)) is the average number of fractal sites separated from the origin by a distance bracketed by r m and 
r m + Ar. This is essentially the starting relationship used (implicitly) by Havlin et al. || to find that, for large N, 

S N (t) ~^/ 2 (lnAO rf// " (4) 

in the non-trivial time regime (or regime II) [Q. Here u — d w /(d w — 1), d w is the anomalous diffusion exponent, df is 
the fractal dimension of the substrate, and d s = 2dj/d w is the spectral dimension. However, the hypothesis leading 
to Eq. (||) is in general false, as we will explicitly show in Sec. [n]by means of numerical simulations of r t (r mj i) for the 
two-dimensional percolation cluster at criticality. Indeed, it is known |)| that the fluctuations of the probability density 
P(r, t) of random walks (also called the propagator or Green's function), which is a statistical quantity closely related 
to the survival probability, exhibits a broad logarithmic distribution for some random fractals such as percolation 
clusters and self-avoiding walks. Bunde et al. Jj| have found that the quantity {P(r, t) q ) exhibits multifractal scaling, 
(P(r, t) q ) ~ (P(r, t)) T ( q \ where r(q) ~ q 1 and 7 = (d e w — l)/(d w — 1) and d l w is the chemical random walk dimension. 
This behavior is a consequence of the large fluctuations of P(r, t) for fixed r and t from a given aggregate to another. 
Nevertheless, these authors have also shown that the distribution of the propagator in the chemical £ space, P(£,t), 
is narrow and, consequently, (P(£,t) q ) ~ (P(£,t)) q . The chemical distance £, the length of the shortest path between 
two sites along lattice bonds, is a more natural measure than the Euclidean distance in disordered systems. It is, for 
example, the distance used in the calculation of optimum paths in cities. 

Let £ m ^ label the i-th site out of those n(m) that are placed at a chemical distance £ from a given origin with 
l m < I < £ m +i, £m = mAl and A£ small (say, of the order of the lattice spacing), and let T t {£ m ,i) be the survival 
probability in the chemical space defined as the probability that site £ m ^ has not been visited by time t by a single 
random walker starting from the origin. Then we can rewrite Eq. (^|) in the chemical £ space as 

00 ln(m) \ 

M*) = E(E {i-irMCOf}) ■ (5) 

m=0 \ i=l I 

One may expect that the distribution of T t (£m,i) for fixed £ m and t is as narrow as the distribution of the propagator 



in the chemical space. (We support this conjecture in Sec. Ill by means of numerical simulations of T t (£ m .i) in 
two-dimensional percolation clusters at criticality.) In this case 

T t {l m<i )^{T t {l mti ))=Y t {l m ) (6) 

for all possible lattice realizations so that ([Tt(£ m ,i)] N ) — (^t(£m,i)) N and, therefore, SW(t) can be approximated by 

00 

S N (t)=J2{l-lr t (£ m )] N }(n(m)) , (7) 

m=0 

where (n(m)) is the average number of fractal sites separated from the origin by a chemical distance with value 
between £ m and £ m + A£. From this formula and following the procedure outlined in the preceding paper pj, in Sec. || 
we will arrive at an expression for S^it) for the non-trivial time regime whose leading asymptotic behavior coincides, 
apart from the value of the prefactor, with the recent proposal, based on a scaling approach, of Drager and Klafter 

m 

S N (t) ~ t d "/ 2 (lnN) d ^ v (8) 

with v = d l w /{d l w — 1). Equation (||) differs from the relationship proposed by Havlin et al. ||, Eq. (|4[), for those 
cases such as that considered in this present paper where <i m i n 7^ 1. Both Havlin et al. and Drager and Klafter 
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supported their conjectures by means of data collapsing plots of computer simulation results obtained for two- and 
three-dimensional percolation aggregates, respectively. We here want also draw attention to the risk involved in this 
method of analysis when the influence of the corrective terms is not properly considered since these terms have a large 
influence jll[] on the final value of Sjv(f). 

The paper is organized as follows. In Sec. |TJ we present the asymptotic evaluation of Sjv'(i) for stochastic fractal 
lattices as a translation to the chemical space of the procedure implemented in the Euclidean and deterministic fractal 
cases . A less rigorous but fairly simple and instructive method for obtaining the main asymptotic term and 



estimating the corrective terms of Sjv(i) for large N is also presented. In Sec. Ill we report simulation results for the 
survival probability of a random walker on a two-dimensional incipient percolation aggregate when a trap is placed 
at a site at a fixed chemical distance I or Euclidean distance r. We find that the distribution is narrow [broad] if the 
traps are located at a fixed chemical [Euclidean] distance. The parameters governing (i) the asymptotic behavior of 
T t (£) (c, v, [a and A), (ii) how the fractal volume grows (Vq and dg), and (iii) how fast a single walker diffuses (Dg 
and d l w ) are estimated in this section. In Sec. |lV| we compare the zeroth- and first-order asymptotic expansion for 
£jv(i) with simulation results obtained for the two-dimensional incipient percolation aggregate. We also criticize the 
reliability of typical collapsing plots for the determination of the dominant trend of SN(t). We conclude with some 
remarks in Sec. 

II. TERRITORY COVERED BY TV RANDOM WALKERS ON A STOCHASTIC FRACTAL SUBSTRATE 

In this section we will translate the results of the previous paper H to chemical language. Reasons for this procedure 
have already been given in Sec. |. We start by replacing Eq. (?]) by its continuum approximation 

S N (t)^J [l-[T t {l)f}d l V$l d *- 1 dl 

« VfaQDtf'W 1 *' J™ {l - {TtiOf}^' 1 ^, (9) 

where dV(£) = Vq dg £ dt ~ 1 d£ is the average number of fractal sites placed at a chemical distance between £ and I + d£, 
and £ = £/(y/2Dgt 1 / d ™). Here Dg is the diffusion constant defined by the Einstein relation, 

L 2 = 2Dgt 2/d ™ , (10) 

where L 2 = (£ 2 ) is the mean-square chemical distance traveled by a single random walker by time t (t large). Next, 
we assume that the asymptotic dependence of T t (£) for £ 3> 1 is given by 



T t (£) « 1 - AC^e-^ (\ + h n r nv j , (11) 

with v = d^j/id^ — 1). This functional form holds (at least in its first terms) on Euclidean lattices BM O , H and 
agrees with the expression conjectured in Refs. [0J^] for fractal substrates with d m i n = 1 such as the Sierpinski gaskets. 
The dominant asymptotic behavior of the propagator in chemical space HHH, P(l,t) ~ exp(— eg"), also coincides 
with the assumed dominant exponential decay of the mortality function 1 — T t (£) in Eq. (pd|). It is known that both 
the propagator and the mortality function share the same asymptotic behavior for Euclidean lattices and for the 
Sierpinski lattice |l4j and we can expect that this coincidence also is the case for stochastic fractals (we will check 



this supposition in Sec. III). The rest of the analysis is identical with that carried out for Euclidean and deterministic 
fractal lattices save for the change of the Euclidean parameters to their chemical analogues. The result for S^(t) 
is, consequently, 

/ , oc n \ 

s N (t) w s N (t) i - ^ Yl lir " N 12 s ™ ) ln '" ln N ) ( 12 ) 



V 

n—1 m— 



with 



S N (t) = V e (2D e ) d ^ 2 t d '/< f — ) (13) 



dg/v 
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and 



*P = A. 



,(2) 



0T 



(2) 
S l = - 

( 2 ) 



/i~ + (f3 — 1)/xcj 



(14) 
(15) 

(16) 

(17) 
(18) 



Here cj = 7 + ln^4 + /ilnc, 7~ 0.577215 is the Euler constant, and — dg/v. The dependence on t and N of the 
main term of Sjv(i) as given by Eq. (13), i.e., ~ i ds / 2 (ln N) dt / V , coincides with the prediction of Ref. fUM. 



A simpler way to estimate the territory covered 



We finish this section by showing how to find the full main term of Eq. (|1J) and even predict the form of the 
corrective terms by only resorting to extremely simple arguments. The ideas of the following reasoning were already 
used in Ref. jl5] ]. The crucial point in our argument is that, for a fixed time t, 1 — [r t (^)] Ar approaches a unit step 
function 0(£ — £ x ) when N — ► oo, £ x being the step's width (see Fig. |l|). The reason for this behavior is clear: For 
large N, [T t (£)] N is only non-neglible when T t (£) is very close to 1. Obviously this occurs when the root-mean square 
chemical distance L{t) traveled by the single random walker by time t is small compared with £, i.e., when £ = £/L(t) 
is large. This in turn implies that in the evaluation of Sjs;{t) only the behavior of T t (£) for large £ is relevant. Then, 
as 1 — \Pt(£)] N approaches a step function of width £ x , the integration of Eq. (||) yields 

S N (t)KVSt**, (19) 

i.e., the territory covered is just the volume of a chemical hypersphere of radius £ x ■ Defining the width £ x of the step 
function as the distance at which 1— takes the intermediate value 1/2 (any other value between and 1 would 

also be valid as £ x is not very sensitive to this value when N 3> 1), and assuming that 1 — w A(;~^ v exp(— c£") 
for large f, we deduce that 1/2 w NA£~^ V exp(— c£x), with £ x = £ X /L so that 

c£ v x w In N - fiv In £ x + In 2 A. (20) 

The term In iV is dominant on the right-hand side of Eq. ( |20| ) for large N, so that a first-order solution of this equation 
is 



c? x «lniV, 



(21) 



i.e., £ v x w L v \yi{N)/c. Hence Eq. ((TJ) yields 



S N (t) « V^£ d| 



In AT 



dpfv 



(22) 



which is in full agreement with the main term of Eq. (|12j) when the Einstein relation, Eq. (|10|), is considered. 
Inserting the above first-order solution for £ x into the right-hand side of Eq. (po|), we get the improved solution 
c£ v x rs In AT — In In AT + In Ac^ + In 2, so that Eq. (O) becomes 



S N (t) 



In AT 



i 



d e -/^ In In AT + lnAc^ + In 2 



In AT 



(23) 



This expression is strikingly close to the first-order approximation of Eq. ( |12| ) , the only difference being that the term 
In 2 = 0.693 • • • in Eq. ( p3| ) plays the role of the Euler constant 7 ~ 0.577215 in Eq. ([l2|). Finally, notice that this 
simple method is not limited to disordered media but that it is also valid for estimating S , at(^) for the non-disordered 
substrates (Euclidean and deterministic fractal media) considered in Refs. 10- 
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FIG. 1. Function 1 - [r t (£)] N versus £ = £/L(t) for (from left to right) TV = 1,10,100,1000 and 10000 where 
Tt(£) = 1 - £~ M " exp(-£"), c = 1.0, u = 1.7, fi = (solid line) and /j = 0.8 (dashed line). We have not plotted the un- 
physical values that appear in the case with n = 0.8 when £ goes to zero. Notice the large influence of the subdominant power 
term on the value of Vt(£)] N which will be reflected in the value of Sjv(i)- The triangles mark the value of £ x for (from 

left to right) N = 10, 100, 1000 and 10000 obtained from Eq. @. 



III. SURVIVAL PROBABILITY, FRACTAL VOLUME AND DIFFUSION IN PERCOLATION 

AGGREGATES 

We have carried out simulations for the number of distinct sites visited by N independent random walkers on a 
typical stochastic fractal: the percolation aggregate embedded in two dimensions. In our simulations every random 
walker makes a jump from a site to one of its nearest neighbors placed at one unit distance in each unit time. The 
incipient percolation aggregates were constructed by the standard Leath method |^|,|l6| on a square lattice with side 
400. In the Leath method, a seed is placed on the site in the center of this box and the cluster is generated by epidemic 
spreading to their nearest neighbors with an infecting probability p c = 0.5927460 || corresponding to site percolation 
in the square lattice. At every step of the generation process a new chemical shell is added to the previous shell by 
the occupation with probability p c of its empty nearest neighbors. The process continues until we reach a shell whose 
sites do not infect any of their empty neighbors; in this case the aggregate so generated is rejected. If the cluster 
grown in this way spans the box from a side to the other, the cluster is accepted as a good representation of a portion 
of an infinite percolation aggregate. Our simulations were carried out over 2000 aggregates generated in this way. 

In order to compare the simulation results for Sn (t) with the predictions of our theoretical approach, Eq. ( |l2| ) , we 
must check that the survival probability or, equivalently, the mortality function, h(£,t) — 1 — T t (£), really behaves 
in the form conjectured in Eq. (pd|). Moreover, we must confirm first that, for a given chemical distance I, the 
distribution of h(£,t) over different realizations of the incipient percolation cluster is narrow because our theoretical 
analysis [cf. Eq. (Q)] relies on this assumption [cf. Eq. (g)]. The numerical evaluation of this quantity as well as 
the propagator P(i, t) (i.e., the probability of finding a single random walker at site i at time t) is performed by the 
Chapman-Kolmogorov method (also called the exact enumeration method HQ). Initially the propagator takes the 
value 1 at the origin site, P(O,0) = 1, and at any other site. This density evolves at every time step by updating 
its value at every site in the form described by the master equation 



where b is the maximum coordination number (6 = 4 for the two-dimensional percolation aggregate) and M is the 
number of neighbors of site i that belong to the aggregate. This evolution equation is valid for the so-called blind 
"ants" ||[l7| because the random walker "attempts" a jump at time t + 1 to a possible nearest neighbor (selected 
randomly) of the site it occupied at time t. If the selected site does not belong to the cluster, the random walker 
stays at the same site and no jump takes place. The "blindness" of the random walkers is taken into account by 
the second term of the right-hand side of Eq. ([24J) . The trap is simulated by a special site belonging to the cluster 




neighbor— 1 



(24) 
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FIG. 2. Plot of the histogram Af(\nh) versus 
values are r = 30, I = 80 and t — 1000. 



\nh\ for fixed r and t (dotted line) and for fixed I and t (solid line). The 



that absorbs all the probability density that enters it without giving back any probability to its neighbors. In the 
simulation of the mortality function, we located a trap at a chemical distance £ = 30 in each of the 2000 percolating 
clusters. We repeated the experiment for traps located at a fixed Euclidean distance, r = 80. The resulting histogram 
for t = 1000 is shown in Fig. ^ (to be compared with the histogram of the propagator shown in Fig. 4 of Ref. Q). 
One observes that the distribution corresponding to fixed £ is very narrow whereas the Euclidean version is broad and 
exhibits a long tail. The striking contrast between the two histograms in Fig. [| can be understood qualitatively. A 
fixed chemical distance I between the origin and the trap means that there exists at least a minimum path connecting 
those sites whose length is t. Thus, the minimum time taken by a random walker to arrive at the trap site is t = £ 
independently of the lattice structure. On the other hand, a fixed Euclidean distance r could correspond to many 
different chemical distances from one cluster to another. This is a consequence of the stochastic lacunarity of the 
fractal aggregate. A large hole between the origin and the trap implies a large minimum time to travel around the 
border of that hole and, obviously, a small mortality in comparison with another cluster where no holes hinder the 
diffusion. 

Figure |^ shows the chemical mean-square displacement 

L 2 = (e 2 ) = /j2e 2 p(£,t)\ (25) 

\ sites / 

as a function of time. The propagator in the chemical space P(£,t) is obtained by summing P(i, t) over all cluster 
sites i on the chemical shell situated at distance £ from the origin. The result is compatible with the Einstein relation 
Eq. © with 2D e = 1.20 ± 0.1 and S w = 2.40 ± 0.05. This value for d e w coincides with that obtained in Ref. § and 
is in agreement with the value reported in Refs. [p|l^]. 

In Fig. [|we plot ln(— \nh(£, t)) versus £ = £/t 1 / d ™ with £ = 80 and, according to the previous discussion, d l w = 2.40. 
If the conjecture in Eq. ( |ll"| ) is right, we can take h{£, t) ~ exp(— c^ v ) as a first approximation, and hence should observe 
the linear behavior ln(— bxh(£, £)) ~ lnc + v£ with £ = \fTD[^ and c = c/(2Di) v ^ 2 . Certainly the plot seems linear 
except for a portion in the range £ > 2.2. This is a finite size effect (already analyzed in the case of the two-dimensional 
Sierpinski gasket in Ref. p4|) associated with the existence of a minimum arrival time corresponding to a random 
walker who travels "ballistically" along a chemical path from the origin to the trap, which in turn implies a maximum 
available value of £ in the simulations (in our simulations this maximum value is SO/SO 1 /' 4 ™ ~ 12.9). A reliable 
interval for numerical fits should exclude this very short time regime. A linear fit in the interval 1.6 < ln£ < 2.17, 
corresponding to 200 <t< 1000, gives the values c = 1.2 ±0.1, i.e., c = 1.3 ±0.1, and v = 1.6 ±0.05. The dashed line 
in Fig. ^ corresponds to these values. The good agreement with numerical values in the above interval seems to assure 
the correctness of the approximation h(£, t) ~ h & (£, t) = exp(— c£, v ) with the values of c and v given above. However, 
the solid line in Fig. [| is a challenge to this interpretation: one sees that the function hb(£, t) — exp(— c£ u ) 
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FIG. 3. Plotofln^ 2 ) versus hit, (I 2 ) = L 2 being the chemical mean-square displacement of a single random walker calculated 
over 2000 two-dimensional percolation clusters at criticality. The line represents the function L 2 = 2Di t 2 ^ dul with 2Di = 1.2 
and d e w = 2.4. 

with v = 1.70, c = 0.9 (i.e. c = 1.05), /i = 0.8 and A = 1.1 is as good as ha,(£,t). Indeed, hb(£,t) is more consistent 
from a theoretical point of view than h a (£,t) because the expected theoretical value of v corresponding to d w — 2.40 
is v = d^/ (d l w — 1) = 1.71, which is in better agreement with the exponent v = 1.7 of h a (£, t) than with the exponent 
v = 1.6 of hb(£,t). Finally, it should be noticed that the values c = 1.05, v = 1.70 are also in agreement with the 
corresponding parameter values of the propagator ||, thus supporting the guess made in Sec. |l| [see below Eq. Q)] 
that the dominant exponential term of the propagator and of the mortality function are the same. This leads us to 
consider that the set of parameters c = 0.9, v = 1.7, fj, = 0.8 is more reliable than c = 1.2, v — 1.6, /i — 0. Obviously, 
further intensive (and extremely time consuming) computer simulations for the mortality function would be required 
in order to reliably determine the values of the parameters that appear in Eq. (jll]) and in the asymptotic corrections 
of S N (t). 

We have also evaluated numerically the fractal volume in terms of the chemical distance V(£), i.e., the number of 
lattice sites inside a circumference (in chemical space) of radius £ . The results are shown in Fig. |5|. A good fit to 
dV{£) = d i V,$£ d < ! ~ x d£ is found with V e = 1.1 ± 0.2 and d £ = 1.65 ± 0.05. Taking into account that d f = 91/48, we 
deduce that d m ; n = df/de — 1.15 ± 0.05, which agrees with previous estimates H[l]. 

IV. SIMULATION RESULTS: TERRITORY COVERED BY TV RANDOM WALKERS ON THE 

PERCOLATION AGGREGATE 

As discussed in Sec. |, the average implicit in the evaluation of S^it) is double: first, we take an average over 
experiments performed on the same aggregate and then a second average over different percolation clusters. The N 
random walkers are always placed initially upon the site in the center of the square box. In our simulations we have 
performed an average over 100 runs per cluster and a second average over 2000 percolation clusters in order to achieve 
good statistics. The maximum time considered was t = 1000. 

According to Eqs. ( [l2| ) and (|l3|), the quotient S n (t) / (In N) 7 with 7 = di/v is only a function of t. In Fig. || the 
logarithm of that quotient is plotted versus lnt for several values of N. The data collapse and the slope close to 
0.66 seems to support Eq. (|^) with 7 = di/v = 0.97, which is in agreement with similar recent results for the three- 
dimensional percolation aggregate Jipj . The collapse is, however, slightly poorer when the exponent 7 = df/u = 1.24 
(df = 91/48 and d w = 2.87 HU) proposed by Havlin et al. [|| [see Eq. (|4|)] is used, as Fig. ^ shows. So, one migth 
be led to the conclusion that the correct value of 7 as defined above is di/v. But in this analysis there was no 
consideration of the relatively large logarithmic corrections predicted by the asymptotic analysis presented in Section 
||, so that the reliability of the above conclusion is seriously affected by this omission. 

To illustrate this point, let us now carry out the same kind of analysis with the simulation results of SN(t) when 
the substrate is a three-dimensional Euclidean lattice. For this case it is well known ^,0] that S'iv(i) is given by an 
asymptotic expression with the form of Eq. (|l2|) in which the logarithmic corrective terms are very important even 
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FIG. 4. Plot of ln( — In h(£, t)) versus ln£ averaged over 2000 two-dimensional incipient percolation clusters. The trap was 
always placed at a site at distance £ — 80 from the origin. The lines represent the functions h(£,t) = exp(— c£ v ) where 

{A = 0.9, /i = 0.8, c = 0.9, v = 1.7} (solid line), {A = 1, fj, = 0, c = 1.2, v = 1.6} (dashed line), and {A = 1, ^ = 0, c = 0.9, 
w = 1.7} (dotted line). 




In ( 



FIG. 5. Plot of \ndV(£) versus \n£ where dV(£) = V(£ + 1) — V(£) is the average number of sites in the chemical shell at 
distance £. The line represents the function Vodil dl ~ 1 with Vq = 1.1 and = 1.65. 
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FIG. 6. Plot of ln(Siv(t)/(lniV) 7 ) versus \nt for N = 2 2 (diamond), 2 5 (down triangle), 2 s (up triangle), 2 11 (circle), and 
2 13 (square) in the two-dimesional percolation aggregate with 7 = de/v ~ 0.97 and 7 = dj/u ~ 1.24. The values corresponding 
to 7 = d«/« ~ 0.97 have been shifted up by 3/2. The line has a slope equal to d s /2 ~ 0.66 



for very large values of N. Indeed, the main asymptotic term leads to very poor predictions for Sjv(t), whereas the 
second-order approximation (n = 2) gives excellent agreement with numerical simulation results. The exponent 7 of 
the main logarithmic term in N and the time exponent dt/d l w — d s /2 are equal to 3/2. We have plotted in Fig. the 
quotient S/v(£)/ In 7 N versus lni for several values of N taking into account that the rigorous value of 7 is 3/2. We 
see that the collapse is far from being perfect because the logarithmic corrections have been ignored. Nevertheless, an 
effective (but incorrect!) value of 7 = 2.75 yields a much better data collapse and a slope close to the theoretical value 
d s /2 — 1.5. We thus conclude that analysis of data collapse plots based on the form of the main term of quantities 
such as S]y(t) (which typically exhibit large corrective terms) should be performed with caution. The values of the 
exponents estimated in this way are untrustworthy because the existence of logarithmic corrections to the main term 
cannot simply be ignored. The value of 7 = 2.75 obtained before is then only an effective way of including all these 
corrective terms together but the true expression involves a main term of the form (ilniV) 3 / 2 times a series similar 
to that given in Eq. (p2J). These considerations should prevent us from drawing hasty conclusions from a simple view 
of plots such as Figs, gfand [?|. 

Finally, in Fig. we show the dependence of Sjv(i) on N and compare simulation results with the zeroth- and 
first-order asymptotic prediction given by Eq. (fl2"|). When the parameter set c = 0.9, v — 1.7, [i — 0.8,^4 = 1 (see 
Sec. Ill) is used, we get results with a very familiar aspect as they are quite similar (although, perhaps the first- 
order approximation is too good) to that already found for Euclidean ||J| and Sierpinski lattices jj] . This is indeed 
encouraging. However, when the parameters c = 1.3 and v = 1.6 are used, we obtain a surprising and strikingly 
accurate zeroth-order approximation. At this point, we again suspect that this last set of parameters are only 
effective parameters that include the influence of the true logarithmic corrective terms in the range of N simulated. 
Hence, Fig. ^ illustrates again, but from a different perspective, how the omission of important corrective terms could 
lead to finding effective parameters that, although providing excellent approximations in the (relatively short) range 
under consideration, are really erroneous. 



V. SUMMARY 



In this paper, the average fractal territory covered up to time t by N independent random walkers all start- 
ing from the same origin on stochastic fractal lattices is calculated in terms of an asymptotic series expansion, 
E~=oEm=o s «™( lnA d " A '~ n ( lnlnA m [see Eq. (||)], which is formally identical to those obtained for Euclidean 
and deterministic fractal lattices. Equation ( p"2| ) is obtained by assuming that (i) the average fractal volume inside 
a "hypersphere" of chemical radius r grows as Vtfr de , (ii) the distribution of the the short-time survival probability 
of a single random walker in the presence of a trap is narrow, so that Eq. (|^) holds, and (iii) this short-time sur- 
vival probability is asymptotically given by Eq. (|Tl"|). We performed numerical simulations for the two-dimensional 
percolation aggregate at criticality which support the validity of the above assumptions. The zeroth- and first-order 



1 

2.5 3.0 3.5 4.0 4.5 5.0 5.5 
In t 

FIG. 7. Plot of ]n(S N (t)/Qn iV) 7 ) versus hit with N = 2 8 (diamond), 2 10 (down triangle), 2 12 (up triangle), 2 14 (circle), 2 
(square) for the three-dimensional Euclidean lattice with 7 = 3/2 and 7 = 2.75. The line has a slope equal to 1.4. 
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FIG. 8. Plot of SN(t)/t de ^ d ™ versus In N in the two-dimensional incipient percolation aggregate for N = 2°, 2 1 , . . . , 2 13 . The 
circles [squares] are the simulation results for t — 1000 [t = 500] averaging over 2000 aggregate realizations. The dashed [dotted] 
line is the zeroth-order theoretical prediction with c = 1.05 and v — 1.7 [c = 1.3 and v = 1.60] and the solid line is the first-order 
approximation with c = 1.05, v = 1.7, // = 0.8 and A = 1. 
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theoretical asymptotic expression of Sjv(i) were calculated explicitly and, for the two-dimensional percolation aggre- 
gate, they compared reasonably well with numerical simulation results. The agreement is similar to that found for 
non-disordered media. 

In our procedure, the use of the chemical distance turns out to be fundamental because the distribution of the 
short-time survival probability in the chemical space is so narrow that we can safely replace the power N of the mean 
value of the survival probability by the mean value of the power N of the survival probability. (On the contrary, this 
does not hold at all when Euclidean distances are used.) This allowed us to easily translate the theoretical results 
previously derived for Euclidean and deterministic fractals to disordered media. 
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